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Abstract 

The neutrino signal from SN 1987A provides an excellent opportunity to constrain 
ly-j physical theories for matter at extreme conditions and properties of particles that are 

<3\ produced in supernova (SN) cores. Phase transitions in the supranuclear equation 

of state (EOS) may change the cooling history of the collapsed stellar core and may 
even lead to a collapse of the protoneutron star to a black hole on a time scale of 
several seconds. The implications of such a scenario for the characteristics of the 
neutrino emission are tested. Moreover, the consequences of a reduction of the neu- 
trino opacities by a possible suppression of the axial-vector interaction in a dense 
■^j- medium (Raffelt and Seckel 1994) are investigated. Since the strongest constraints 

C*~) from SN 1987A for the considered theories rely on the signal duration marked by the 

last three events in Kamiokande, we also discuss attempts to explain these events dif- 
ferent from the standard protoneutron star cooling scenario. In particular, attention 
is focussed on the scenario of a rapid cooling of the protoneutron star within the time 
interval of the first 8 Kamiokande events or the 1MB signal which is only about 2-5 s. 
Such a short neutrino emission could be the consequence of low neutrino opacities or 
convectively enhanced neutrino transport. Since the driving force of neutrino heat- 
O ing around the protoneutron star diminishes, the pressure support of the overlying 

£j material decreases and some matter might fall back to the neutron star. However, 

analytical estimates and numerical simulations indicate that the neutrino emission 
associated with the accretion of this material is hardly energetic enough to account 
for the last bunch of neutrinos at Kamiokande. 
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1. Introduction 

Massive stars (M > 8M©) develop iron cores at their centers during a last hydrostatic 
burning phase that uses Si as fuel. When the mass of this core gets close to the Chan- 
drasekhar mass, electron captures onto protons and nuclei and, for very massive stars, 
also photodisintegrations of nuclei speed up and reduce the pressure. Therefore the 
core becomes unstable against gravitational collapse. Within a fraction of a second 
the collapse proceeds to nuclear matter densities, where it is finally stopped due to the 
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stiffening of the equation of state by repulsive nuclear forces and the pressure of non- 
relativistic (degenerate) nucleons. At this moment a strong shock wave is launched and 
starts to propagate outward in mass and radius. However, this shock is unable to cause 
the disruption of the progenitor star directly. Its propagation into the stellar mantle 
and envelope, which then leads to the spectacular outburst of light in a supernova ex- 
plosion, is driven by the deposition of energy by neutrinos in the material behind the 
shock. These neutrinos are emitted in very large numbers from the hot protoneutron 
star that is formed at the center of the explosion and that releases its gravitational 
binding energy (several 10 53 erg) in neutrinos. 

On Feb. 23, 1987, at 7:35:35 UT(±lmin) a bunch of 11 neutrinos within a 12.4 s 
time interval was recorded by the Kamiokande II collaboration (Hirata et al. 1987) and 
another 8 neutrinos were counted in the 1MB detector between 7:35:41.374 UT and 
7:35:46.956 UT (Bionta et al. 1987). The Russian experiment at Baksan reported 5 
events (Alexeyev et al. 1987) but has a rather large background rate and could uncover 
the connection with the other two detections only later due to problems with the time 
measurement. The neutrino events coincided in time with SN 1987A in the Large Magel- 
lanic Cloud, which was first observed only a few hours later. The overall characteristics 
of these neutrino detections like the event energies, the number of events, and the burst 
duration, agree well with theoretical predictions for the neutrino emission from stellar 
core collapse events. Both facts leave nearly no doubt that the neutrinos came from 
SN 1987A and thus mean the first detection of neutrinos from an extragalactic source. 

This historical detection was considerably useful to constrain properties of neutrinos 
and of other particles that are possibly produced at the conditions in hot supernova 
cores. Moreover, the neutrino burst means a strong confirmation of principle aspects of 
supernova models and of the theoretical picture of the very complex processes that occur 
during supernova explosions. Certainly, the few neutrino events provide rather poor 
statistics and (therefore?) the recorded data show a number of unsatisfactory features, 
i.e., the weak compatibility of the Kamiokande and 1MB signals, the clear peaking of 
the angular distributions away from SN 1987A, although isotropy was expected, and 
the 7 s gap between the eighth and the ninth event in Kamiokande. Nevertheless, the 
SN 1987A neutrinos have to be considered as a serious test for new physical models of 
matter at high densities or particle processes in a dense medium both of which may 
affect the evolution of supernovae and of new-born neutron stars. A nearby Galactic 
supernova will yield a much larger number of events. Hopefully, a more precise time 
measurement in different experiments and a much better temporal resolution will give 
us more detailed information about the spectral and temporal signature of the neutrino 
emission. 

Although the neutrino signal yields the most direct constraint on the physical 
processes deep inside the exploding star, it is not the only one. Further constraints can 
be deduced from observed properties of neutron stars (masses, kick velocities, cooling 
behavior), from supernova nucleosynthesis, which was recognized during the past few 
years to depend sensitively on the character of the neutrino emission from the nascent 
neutron star, and from the explosion energy of a supernova, which is determined by the 
neutrino-matter interaction near the protoneutron star. Moreover, also the structure 
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of the photon flux from a supernova provides information about the processes during 
the onset of the explosion. SN 1987A showed that the light curve is strongly affected 
by large-scale mixing which reaches deep into the star and which occurs already during 
the early moments of the explosion. Radioactive elements that are formed close to 
the center are therefore mixed into the He- and H-layers of the progenitor star with 
very high velocities and with significant anisotropy. This was deduced from infrared 
emission lines of Fe in the spectra of SN 1987A and from X-ray and 7-ray observations 
of SN 1987A and other supernovae. The variety of different observables which are 
influenced or determined by the physical events in the deep interior will make Type-II 
supernovae a very suitable high-energy physics laboratory once a standard model of 
supernova explosions has been developed. 

In the presented paper the main results of two projects are summarized. In the 
first project we (Keil and Janka 1995) investigated the implications of hadronic phase 
transitions at supranuclear densities for the deleptonization and cooling of protoneutron 
stars. These phase transitions are a consequence of the production of new baryonic or 
mesonic states like hyperons, pions or kaons (Sect. 2). In the second project, described in 
Sect. 3, we (Keil, Janka, and Raffelt 1995) studied the consequences of a reduction of the 
neutrino opacities due to a possible suppression of the axial- vector currents in a dense 
medium (Raffelt and Seckel 1994; Janka, Keil, Raffelt, and Seckel, in preparation). In 
Sect. 4 some issues of non-standard models for the neutrino emission from supernovae 
will be discussed with respect to the question whether the 7 s gap followed by the 
last three Kamiokande events might have a physical explanation rather than being 
considered statistical fluctuation. 



2. Hadronic phase transitions at supranuclear densities 

At high densities the Fermi momenta of nucleons can become sufficiently high to allow for 
the production of additional hadronic states besides n and p. The physics at densities 
above nuclear matter density (/9 nuc ~ 2.5 • 10 14 g/cm 3 ) is not well understood and a 
variety of suggestions for the occurrence of new particles can be found in the literature. 
Dependent on the assumptions about the nucleon-nucleon interaction the production of 
pions, associated with the existence of a 7r-condensate, has been considered (see, e.g., 
Takahara et al. 1991, Takatsuka et al. 1992), the formation of a condensate of K~ was 
discussed (Brown et al. 1992, Mayle et al. 1993), and the effects of additional hyperonic 
states (A, S,...) or A-resonances have been studied (e.g. Glendenning 1985, 1989). 

Using finite temperature extensions of two T = EOSs provided to us by N.K. Glen- 
denning (personal communication) we investigated the implications of a hyperonization 
for the neutrino cooling of newly formed neutron stars. EOS A describes a gas con- 
sisting of n, p, e 1 * 1 , and while EOS B contains hyperons and A-resonances that are 
formed at p ^ 2/? nuc via reactions like 

N + N — > N + A + K (1) 

where the K mesons subsequently decay into 7, /i, and v. The occurrence of new 
hadronic degrees of freedom at high densities "softens" the EOS, i.e. flattens the pressure 



3 



1000 F r 




1014 10 15 

p [g/cm 3 ] 

Fig. 1. Pressure as a function of baryon rest mass density, p = n^m^ (m u = 
1.66 • 10~ 24 g), for /j, Ve = and T = 0.5 MeV. 
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Fig. 2. Electron chemical potential as a function of baryon density p for 
p Ve = and T = 0.5 MeV. 
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increase P{p) as a function of density p (Fig. 1) because the nucleon degeneracy is 
reduced and degeneracy energy (or thermal energy), which is relevant for the pressure, 
is consumed in the production of massive particles. By reactions like those in Eq. (1) 
nucleons, in neutronized matter primarily neutrons, are converted into hyperons. This 
reduces the neutron chemical potential p n and leads to a maximum of the electron 
chemical potential p e = p n — p p of about 200 MeV in case of the hyperonic EOS B 
(for p Ve = 0, i.e. deleptonized conditions) (Fig. 2). Although EOS B is a very special 
example and the inherent assumptions about physical parameters and processes may 
be largely different from the preferred choices and the use in other supranuclear EOSs, 
the behavior of P{p) is qualitatively shared with a variety of EOSs (e.g., compare with 
Brown et al. 1992, Takahara et al. 1991, Takatsuka et al. 1992). 




p c [ 1 15 g/ cm 3 ] 

Fig. 3. Baryonic masses of neutron stars versus the central baryon rest mass 
densities p c of the stars. The two displayed cases for each EOS correspond to 
the lepton-rich (Y\ ev = 0.3) and deleptonized {p Ve = 0) situations. In all cases 
T = 0.5 MeV in the whole star was assumed. 

The use of EOS B is also representative in its consequences for neutron star prop- 
erties and the evolution of newly formed neutron stars. The effects of different EOSs 
can be visualized at their M(p c ) relations which relate the neutron star masses to their 
central densities. Figure 3 displays the baryonic masses M&(/9 C ), defined as the total 
number of baryons iV& times the atomic mass unit m u , for the cold (T 0) EOS A and 
EOS B both for the lepton-rich (lepton fraction Y\ ev = Y e -\-Y Ve = 0.3) and deleptonized 
{pv e = 0) cases. The parts of the curves with negative gradients dM&/d/9 c are unstable 
branches, i.e. the corresponding models are unstable against perturbations and cannot 
escape collapse on a dynamical time scale. Differences between the lepton-rich and 
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deleptonized cases start to become significant for p c <; 5 • 10 14 g/cm . EOS A allows 
for neutron stars with larger maximum masses than EOS B. For EOS A the maximum 
baryonic mass is 2.02 Mq (maximum gravitational mass 1.84 Mq) in the lepton-rich case 
and 2.26 Mq (1.96 Mq) for deleptonized configurations. EOS B has a limiting baryon 
mass of 1.94 Mq (1.77 M©) for Fi ep = 0.3 and 1.77 M (1.58 M ) for fi„ e = 0. Although 
the creation of hyperons consumes energy and converts degeneracy or thermal energy 
into particle rest mass energy (thus moderating the pressure increase with density), the 
hyperonization usually is associated with a net energy release and a higher neutrino loss 
of the cooling protoneutron star compared to an equal-mass star with EOS A. This is, 
because the softer EOS leads to a more compact structure and a stronger gravitational 
binding. 

It is particularly interesting to compare the evolution from the lepton-rich state of 
a collapsed stellar core to the deleptonized, final neutron star for EOS A and EOS B. 
The neutronization process, of course, conserves the total baryon number and therefore 
proceeds along a horizontal line in Fig. 3. Starting from a stable state for the lepton-rich 
EOS A such a track always leads to a stable deleptonized configuration, while for EOS B 
the situation can occur that the baryon mass of the collapsed stellar core is higher than 
the maximum stable mass for the deleptonized case. This will inevitably lead to a 
collapse of the cooling and neutronizing remnant to a more compact state, probably to 
a black hole. The collapse does not happen immediately but on a secular time scale of 
several seconds, because it is triggered by the gradual decrease of the lepton content of 
the neutrino-radiating, hot protoneutron star. 

We evolved a set of protoneutron star models from the hot, lepton-rich state shortly 
after core collapse to the final, deleptonized, cold neutron star stage to determine the 
time scale until the collapse occurs and to investigate the corresponding neutrino signal 
(Keil and Janka 1995). The models had different masses and different initial temper- 
atures and were constructed and evolved with either EOS A or EOS B. Model names 
beginning with "S" mean small baryonic mass, "L" denotes large mass (in case of EOS B 
above the maximum stable mass of cold neutron stars), "A" and "B" distinguish be- 
tween use of EOS A or EOS B, and "C" and "H" indicate cool or hot initial models, 
respectively. 

Figure 4 shows the central densities function of time. One can clearly see 

the difference of stable and unstable models. The central densities of the latter increase 
steeply when instability is approached while the stable models develop to finite values 
of p c for late times. All models reveal an initial drop of the density which is associ- 
ated with a phase of increasing temperatures when the protoneutron star deleptonizes 
and degenerate electrons are converted into neutrinos which then downscatter and heat 
the stellar medium. Since the production of hyperons is favored as neutrons become 
more abundant in the star, instability occurs as a consequence of lepton number loss 
by v e emission. At the time of gravitational collapse the neutron stars are still hot. 
Accordingly, the stable models continue to radiate neutrinos with considerable fluxes 
for another ten seconds or even longer (Fig. 5). Since internal energy contributes to 
the gravitational potential and increases the star's gravitational mass, we could not 
find a significant stabilizing influence of thermal pressure once the leptons were gone 
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Fig. 4. Central baryon rest mass densities for the sample of computed pro- 
toneutron star cooling models versus (universal coordinate) time. 
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Fig. 5. Total neutrino luminosities for the sample of computed protoneutron 
star cooling models versus (universal coordinate) time. The luminosities are 
measured at the neutron star surfaces and include contributions from all neu- 
trino types. 
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and hyperonization had softened the supranuclear EOS. The additional mass that could 
typically be supported by thermal pressure was about 0.1-0.2 Mq. Once the unstable 
models start to contract and p c increases, their neutrino luminosities drop appreciably 
because of growing opacities and gravitational redshift. At the moment gravitational in- 
stability sets in and the increase of the central density accelerates, the neutrino emission 
is terminated and no late outburst of neutrino radiation occurs. 

Typically, the collapse happens after less than 10 s of neutrino cooling. Model LBC* 
is an exception to this because both its baryonic mass (1.79 Mq) and its gravitational 
mass (1.67 Mq) were chosen only little above the maximum masses stabilized by EOS B 
for /j, Ve = 0. For all other unstable models the gravitational masses were considerably 
higher above the mass limit of EOS B. Therefore only model LBC* could escape the 
onset of instability for nearly 17 s. Hot, collapsed stellar iron cores, however, have not 
yet emitted much of their gravitational binding energy in neutrinos and are usually 
weakly bound. This means that for "realistic" protoneutron star models the initial 
gravitational mass is not much different from the baryonic mass. The corresponding 
model LBH has a gravitational mass considerably in excess of the maximum stabilized 
mass of EOS B. Model LBC, which also collapses in less than 10 s, has a baryonic mass 
that is about 0.04 Mq above the limit, in contrast to LBC* where it is only 0.02 Mq. 

Since the deleptonization and gravitational instability of our models occurs in a 
time shorter than the signal duration at Kamiokande, we are able to draw the following 
conclusions. If the last 3 Kamiokande events between t = 9.219 s and t = 12.439 s (after 
the gap of about 7 s) are real events and do not have some non-standard explanation 
different from the commonly assumed protoneutron star cooling scenario, then either 
the mass of the neutron star in SN 1987A is below the limiting mass of EOS B (M& < 
1.77 Mq) or a hyperonization and softening of the supranuclear EOS as described by 
EOS B is excluded on grounds of the detected neutrino signal. Of course, this only 
concerns the physical picture and the implications associated with EOS B (which is 
based on an effective relativistic field theory of interacting hadrons; Glendenning 1985) 
but does not mean that one has to discard the possibility of a formation of new hadronic 
states or even of hyperons in general. Also, we found that the collapse to a black hole 
happens "quietly" in the considered scenario of a continuous hyperonization, i.e. it is 
not accompanied by a final, energetic burst of neutrino radiation which might help to 
explain the peculiarities of the Kamiokande data. 

3. Neutrino opacities and the SN 1987A neutrino signal 

Nucleons in the neutron star medium are not isolated particles but interact with each 
other and collide frequently. A number of effects are discussed in the literature, for ex- 
ample many-body corrections like ion-ion correlations (Itoh 1975), nucleon correlations 
(Sawyer 1989) or Coulomb correlations of electrons and electron screening (Horowitz 
and Wehrberger 1991), which may lead to a reduction of the neutrino opacity in a dense 
medium compared to the opacity of a gas of independent nucleons. Recently, Raffelt 
and Seckel (1994) argued that the collisions among nucleons which interact via a spin- 
dependent force cause the nucleon spins to fluctuate fast relative to the other relevant 
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Fig. 7. Expected event numbers and mean energies of captured electron an- 
tineutrinos in the KII and 1MB detectors versus protoneutron star mass. 
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time scales like the scattering time scale of neutrinos which is proportional to the in- 
verse of the energy transfer. If the spin of a given nucleon flips several times during the 
characteristic interaction time, the effective nucleon spin "seen" by the neutrinos will 
be reduced by an averaging effect. Therefore, the axial- vector currents that mediate 
the spin-dependent part of neutrino-nucleon interactions might be nearly completely 
suppressed. Similar to Bremsstrahlung emission of neutrino pairs or axions this effect 
requires the presence of a bystander nucleon in a neutrino-nucleon scattering process, 
i.e. instead of the simple scattering v + N — > v + N Raffelt and Seckel (1994) considered 
the reaction 

v + N + N — > v + N + N. (2) 

Since axial-vector interactions with the axial-vector coupling constant Ca yield the 
dominant contribution to the neutral-current neutrino-nucleon scattering cross section, 

°*n = ^(C$ + 3Cl)-el, (3) 

(and are also the dominant part in charged- current /3-processes) this would imply a 
significant reduction of the neutrino opacity at high densities. In Eq. (3) Cy is the 
coupling constant for the vector current, Gp the Fermi constant, and e v the initial 
neutrino energy. The described effect adds to a possible suppression of the vector and 
axial-vector currents in a bulk nuclear medium, which leads to the reduction of the 
coefficients Cy and Ca relative to their vacuum values (see Raffelt and Seckel 1994). 

In a pilot study we investigated the implications of the conjectured suppression of 
axial- vector currents for the cooling and deleptonization of newly formed neutron stars. 
We simulated the neutrino cooling for a large set of protoneutron star models with dif- 
ferent masses, initial temperatures, EOSs (EOS A and EOS B), and neutrino opacities. 
The neutrino emission of these models was evaluated for the expected neutrino events 
in the Kamiokande II and 1MB detectors. Demanding compatibility of the computed 
neutrino signals with the SN 1987A data allowed us to deduce limits on the strength 
of the suppression. For the purposes of this exploratory study it was sufficient to im- 
plement an opacity modification in a rather schematic way: we multiplied the standard 
axial-vector opacities with a suppression factor and combined the standard opacities 
with these suppressed opacities by weight factors (1 — a) and a, respectively, i.e. we 
replaced C\ in Eq. (3) by 



C'X X < (1 - a) + a 



1 r 10MeV\ 

+ 3-10 13 g/cm 3 / V T J 



(4) 



a = corresponds to the standard case of no suppression, a = 1 yields full suppres- 
sion, and with a < we could also investigate the case of enhanced cross sections (for 
completeness). 

Figure 6 shows the variation of characteristic times as a function of the neutron star 
mass (baryon mass). The standard values for neutrino opacities including fermion phase 
space blocking effects (see Keil and Janka 1995) were used in all models. Filled squares 
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Fig. 8. Characteristic times of neutrino emission and expected signal durations 
in the KII and 1MB detectors versus the parameter a. a = corresponds to 
the standard neutrino opacities (see Fig. 6) and a = 1 means full suppression. 
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Fig. 9. Expected event numbers and mean energies of captured electron an- 
tineutrinos in the KII and 1MB detectors versus the parameter a. a = 
corresponds to the standard neutrino opacities, a = 1 means full suppression, 
and a < belongs to the (hypothetical) case of enhanced opacities. 
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correspond to models computed with EOS B, open squares to those with EOS A. t^i is 
the deleptonization time, t c \ the cooling time, tgo is defined as the time interval after 
which 90% of the total energy emitted in neutrinos have left the star. £kii denotes the 
time interval in which 90% of the predicted events in the Kamiokande II detector occur, 
and tiMB gives the corresponding time for the 1MB experiment. In Fig. 7 the average 
energies of captured v e in both detectors, (£j> e )Kll and (£i> e )iMB, and the expected 
integrated numbers of events iVxn and iViMB are displayed. These predictions have to 
be compared with the measured data of the SN 1987A neutrino signal: t^- 7 ^ 12.4 s, 
t*™ « 5.6s, ( £i , e )Kn ~ 17.4MeV, ( £i , e )iMB ~ 33.9MeV, N$£ = 11, = 8. 

The average event energies in both detectors show remarkable stability, which prob- 
ably reflects a kind of self-adjustment of the position of the neutrino sphere: higher tem- 
peratures in the neutron star (at a certain value of the density) cause higher opacities 
and therefore the neutrino sphere moves outward to a lower density where also the tem- 
perature is lower. In contrast, the characteristic times and event numbers increase with 
the neutron star mass and thus with the gravitational binding energy that is released 
in neutrinos. The flattening of iVxn and iViMB and the slight decrease of trMB towards 
the right ends of the curves indicate the growing importance of gravitational redshift 
for massive neutron stars. This is particularly important for the 1MB detector because 
of its high threshold energy of about 19MeV and its slowly rising detection efficiency 
which make the detector very sensitive to the high-energy tail of the neutrino spectrum. 
Comparing the predicted signals with the SN 1987A data the best agreement is found 
for neutron stars with EOS B and a baryon mass between 1.5 and 1.6 M©. Therefore we 
used our 1.53 Mq model as a reference model for the calculations with varied neutrino 
opacities. 

Figure 8 displays the variation of the times with increasing suppression of the 
neutrino opacities. Filled circles mark models where the axial- vector part in neutral- 
current scatterings on nucleons and in charged- current /3-processes v e + n <->■ p + e~ and 
+ p n + e + were changed in accordance with Eq. (4). Open circles indicate models 
where only neutral-current neutrino-nucleon interactions were modified, and asterisks 
denote the case where the critical density in Eq. (4) was taken as 10 14 g/cm 3 instead 
of 3 X 10 13 g/cm 3 . Figure 9 provides the results of expected event numbers and mean 
energies of captured neutrinos for the whole range of tested suppression parameters a. 
While all characteristic times drop drastically in case of reduced neutrino opacities, the 
event numbers and average event energies show the opposite tendency. Suppression of 
the neutrino opacities moves the neutrino sphere to higher densities and therefore higher 
temperatures. This "hardens" the emitted neutrino spectrum and leads to an increase 
of the predicted number of events. Note that due to the very steep density decline in 
the protoneutron star surface only very little matter is at densities below 3 X 10 13 g/cm 3 
(10 14 g/cm 3 ) and the optical depth of this layer is rather small. 

For a — > 1 both decreasing signal durations and increasing numbers and mean 
energies of the neutrino events lead to a clear discrepancy with the SN 1987A data. 
With a lower or higher protoneutron star mass or a different EOS not all aspects of this 
discrepancy can be compensated simultaneously. If the last 3 events in Kamiokande are 
interpreted as a part of the cooling signal of the nascent neutron star, we therefore con- 
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elude that it is very unlikely that the neutrino opacity of the neutron star in SN 1987A 
was less than about 50% (a <; 0.5) of its standard value. From this limit on a possible 
suppression of axial- vector currents also a limit on the spin fluctuation rate in the nu- 
clear medium can be deduced. The spin fluctuation rate seems to be of the order of the 
thermal frequency ~ T and not as large as suggested by naive perturbative estimates 
which yield values of the order of 20-50 T for typical conditions in protoneutron stars. 

As discussed in detail by Keil, Janka, and Raffelt (1995) and by Janka, Keil, Raf- 
felt, and Seckel (in preparation) the structural similarity of Bremsstrahlung emission 
of axions and neutrino-nucleon scattering in the presence of a bystander nucleon is ex- 
pressed by the fact that the calculation of the rates of both processes involves the same 
spin-density structure function for the nucleons. The form of this structure function de- 
pends on the spin fluctuation rate. Upper limits on the spin fluctuation rate therefore 
translate into limits on a possible suppression of the axion emission rate as expected 
from the Landau-Pomeranchuk-Migdal effect (Raffelt and Seckel 1991). Due to the 
spin fluctuation rate being only of order T the axion emission rate in a dense medium 
appears to be not as much reduced as would be the case for the naive perturbative 
estimate of the spin fluctuation rate. However, it is roughly a factor of 10 lower than 
assumed in previous models for the axion emission from supernova cores. Therefore the 
limits on the axion-nucleon Yukawa coupling are weakened and the bounds on the axion 
mass become about a factor of 10 less stringent, giving m a ^ 10 -2 eV instead of the 
previously claimed limit of m a ^ 10 -3 eV. 

4. The last 3 Kamiokande events — evidence for a non-standard scenario ? 

In the studies summarized in Sects. 2 and 3 the strongest constraints and limits were 
obtained from the duration of the SN 1987A neutrino signal. These constraints relied 
on the signal length in the Kamiokande detector as a measure of the time interval over 
which the major part of the neutrino emission of the nascent neutron star took place. 
Since the signal duration is determined by the last 3 Kamiokande events the strange 
time structure of the Kamiokande signal (as well as the marginal compatibility of the 
1MB and Kamiokande detections) requires a closer investigation of the question whether 
the gap of about 7 s and the following bunch of 3 events within 3.2 s indicate some non- 
standard scenario. If events 9-11 could be convincingly explained without considering 
them as part of the cooling emission of the protoneutron star the arguments used in 
Sects. 2 and 3 would be significantly weakened. 

Suzuki (1989) performed a statistical analysis of the significance of the time struc- 
ture of the Kamiokande signal. Assuming an exponential decay of the event rate in the 
detector he showed that for all choices of the decay time scale and initial event rate the 
probability of a 7 s gap followed by 3 events never gets higher than 1-2%. The maximum 
probabilities were obtained for rather long decay time scales, longer than about 15 s. 
Although being quite small these probabilities are not so small that they are imperative 
for a non-standard scenario, in particular not because 1MB observed two events in the 
7 s gap. It may also well be that a different model of the neutrino emission with a higher 
neutrino flux at early times and a less steep than exponential decline later increases the 
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probability of the Kamiokande signal structure. Certainly, however, considerations of 
scenarios alternative to the simple "standard" protoneutron star cooling model are not 
unreasonable speculations. 

In Sect. 2 we have already mentioned that phase transitions in the supranuclear EOS 
with a slow and gradual formation of new hadronic states, triggered by the progressing 
neutronization of the matter, are not able to increase the probability of the last 3 
Kamiokande events because no late outburst of neutrino radiation occurs, even not in 
the case of a collapse to a black hole after about 10 s of protoneutron star cooling. But 
could accretion onto a black hole that forms at some time between 1.9 s (last event 
before Kamiokande gap) and 9.2 s (first event after gap) account for the bunch of the 3 
late events? The delayed formation of the black hole could happen in a scenario similar 
to the one discussed in Sect. 2. Or could even fallback of matter to the neutron star 
yield sufficient neutrino emission to cause 3 events in the Kamiokande detector? In 
the latter scenario the protoneutron star might release its binding energy within a short 
time of only 2-6 s, the upper end of this time interval being set by events 7 and 8 in 1MB. 
Such a short cooling time could mean a stronger suppression of the neutrino opacities 
than inferred from the analysis in Sect. 3. Due to the disregard of the last bunch of 
Kamiokande events and the correspondingly lower total energy required to explain the 
detections, this would also imply a somewhat smaller (initial) neutron star mass than 
estimated from the whole 12 s signal. 

What are the implications concerning the 1MB signal when the 7 s gap in the 
Kamiokande data is assumed to have a physical rather than statistical cause? Not 
dismissing the last two 1MB events as background, the break in the neutrino detection 
of both experiments is actually only about 3.6 s. This is not much longer than the 
pause between the sixth (at 2.684 s) and seventh event (at 5.010 s) in 1MB. Certainly, 
one has to ask whether it is imperative to search for alternative models of the neutrino 
emission from SN 1987A to explain a gap of only 3.6 s. Yet, the energies of the last 
two 1MB counts (at 5.010 s and 5.582 s) are only slightly above the detector threshold 
(19 ± 5 MeV and 22 ± 5 MeV). The low detection efficiency of 1MB in this energy range 
requires a sizable neutrino flux which should have caused more than 100 neutrino events 
in the Kamiokande experiment that has a much lower threshold of about 7 MeV. The 
absence of a corresponding signal at Kamioka was interpreted as an indication that the 
last 2 1MB events might have occurred due to a statistical fluctuation or their energies 
were estimated systematically lower (Suzuki 1989). Dropping events 7 and 8 of 1MB as 
background would certainly be a radical and rather unsatisfactory step. Keeping the last 
2 1MB counts as real, however, lowers the need for a non-standard emission model and 
moreover will not allow one to improve the weak agreement between the Kamiokande 
and 1MB measurements even with an alternative neutron star cooling model. 

Can a non-zero neutrino mass and the corresponding flight delay of less energetic 
neutrinos account for the gap in the Kamiokande data (and the ~ 2.7 s break in 1MB)? 
For two neutrinos of mass m v and energies e\ and £2 with m^c 2 <C £1, £2 one gets for 
the difference of the flight times 
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which is positive for s\ > £2- D is the distance to the source. Plugging in representative 
numbers for the SN 1987A neutrinos, s\ = 30 MeV and £2 = lOMeV and using D 
50 kpc as the distance to the Large Magellanic Cloud one obtains 



2' 



At 12 = 0.23 s ^ . (6) 

10 eV 2 V ; 

In order to explain a gap of several seconds between a bunch of energetic neutrinos that 
arrive earlier than a bunch of less energetic neutrinos a neutrino mass larger than 10 eV 
has to be assumed, which is excluded for electron neutrinos by current laboratory limits. 

Let us therefore return to the accretion scenario already addressed above. The 
energy of a v e burst that causes a number iVd e t of events v e + p — > e + + n in a water 
Cherenkov detector can be estimated to 

E ~ IT nr', r , 1 ^ (^/50 kpc) 2 

E Va ~ 1.7x10 erg £ {T ^ /5MeY) (MlW2140t) " (?) 

£ is the detector efficiency, averaged over the neutrino spectrum, Mh 2 o is the fiducial 
mass of the detector (2140 1 in case of Kamiokande, 6800 1 for 1MB), and T„ e is the 
spectral temperature of the v e burst, the spectrum being described by a Fermi-Dirac 
distribution with zero degeneracy. The last iVdet = 3 Kamiokande events have T„ e = 
1.6 MeV (see Suzuki 1989). For £ = 1.0 one finds a lower limit on E„ e} E„ e <; 1.6 x 
10 52 erg. In a more careful evaluation Suzuki (1989) estimated E„ e 51^° X 10 52 erg. 
For the bunch of the first 8 events in Kamiokande (T^ e 3.3 MeV) one gets E„ e <; 
2 X 10 52 erg (with £ = 1), which is roughly the same number as inferred from the 8 
1MB neutrinos, Eff B 2.4l2:2 x !0 52 erg (see, e.g., Suzuki 1989, Janka and Hillebrandt 
1989). Assuming equal energies emitted in all 6 kinds of neutrinos the total gravitational 
binding energy released by the collapsed stellar core should be 1-2 X 10 53 erg. This leads 
to a somewhat lower estimate of the (initial) mass of the neutron star in SN 1987A, for 
stiff EOSs between about 1.0 M© and 1.5 M©, for soft EOSs only between 0.8 Mq and 
1.2 M Q . 

In the accretion process onto a compact object with mass M, radius R } and redshift 
e *( B ) the accreted matter emits a fraction / of its gravitational binding energy in 
neutrinos. The mass that has to fall onto the compact object to release an observable 
energy E„ e in electron antineutrinos therefore is 

The inequality neglects gravitational redshift and assumes an efficiency of conversion of 
gravitational binding energy into neutrinos of / = 1. The factor f„ e is the fraction of the 
released energy that is radiated as v ei optimistically chosen to be f„ e =1/2 (assuming 
that the energy is radiated only in v e and v e and equally shared between both of them). 

The mass that potentially falls back to the neutron star or black hole at the center 
of the supernova can roughly be estimated analytically, because at times later than 
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about 2 s after shock formation a number of simplifying assumptions can be applied. 
From the virial theorem (Cooperstein et al. 1984, Bethe 1993a) the pressure behind the 
shock front can be estimated as 

^ = t% • (9) 



E sn is the supernova explosion energy, R s the shock radius. The density outside the 
shock is given from the progenitor star by (Bethe 1990, Fig. 22) 

po(r) = Hr~ 3 (10) 

with H sa 6 X 10 31 g for an enclosed mass of 1.8 M Q in the center of the M sa 18-20 M Q 
progenitor of SN 1987A. From the Rankine-Hugoniot jump conditions in the shock one 
gets for the shock velocity 



Us 



p xl/2 , E N l/2 

1 s \ / ^sn 



Po(Rs)J 



1.2 • 10 9 cm/s (11) 



where the postshock pressure and density were assumed to be much larger than the 
preshock values, P s ^> Po(R s ) and p s ^> po(R s ) } and where the expression in the middle 
was obtained by using Eqs. (9) and (10). The numerical value gives the result for 
E sn = 10 51 erg. Note that u s is independent of R s which is not true in the start-up phase 
of the shock but is well fulfilled at later times when the shock is far away from the neutron 
star and pressure gradients as well as gravity have become small (see Bethe 1993b for 
an analytical discussion of this point). In that phase the kinetic energy (per unit mass) 
of the material dominates its gravitational binding energy and further acceleration of 
the matter behind the shock plays no role any more. The explosion energy of the 
supernova E sn is then mostly kinetic energy. In this situation the separation radius r cut 
of supernova ejecta and fallback material can be estimated by equating kinetic energy 
and gravitational binding energy, 

1 9/ , GM , s 
-u 2 r cut = • 12 

2 r cut 

Since according to Eq. (11) the shock velocity is constant, the velocity of the mate- 
rial at a point r should be ~ r/t where t is the time after stellar core bounce and shock 
formation, thus 

/ \ '"cut /-, r>\ 

u(r cnt ) ~ - — r cut — . (13) 

'-cut -Ks 

For the same reason the shock position at the moment of separation, t cut = to + r cut /c s , 
is approximately given by 

Rs u s t cut = u s (t + r cut /c s ) . (14) 

to is the time (after bounce) when the neutron star collapses to a black hole or when 
the neutrino luminosities decrease and the driving force by neutrino heating diminishes. 
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''cut / c s gives a lower limit for the time in which the vanishing of the pressure support at 
the center is communicated outwards and arrives at the radius r cut of the future mass 
cut. This information propagates outwards with a velocity that is somewhat smaller 
(delay by the inertia of the matter) than the sound speed c s in the region behind the 
shock. Using Eq. (13) in Eq. (12) yields 



cut 



2GM (15) 



R s J u 2 s R s 

Assuming constant density in the region enclosed by the shock, which is consistent 
with the assumption of a negligible pressure gradient, Eq. (15) gives an upper limit on 
the fraction of the shocked stellar material that falls back to the neutron star or black 
hole. The fallback mass is thus optimistically estimated as 

AM acc « ^j^J -(M s -M) (16) 

when M s is the mass enclosed by R S} i.e., the mass of the progenitor star which has been 
passed by the supernova shock by the time t cut . If the matter of the progenitor star 
ahead of the shock has not fallen inwards significantly since the onset of core collapse 
(which is well fulfilled at t <; 2 s after bounce because the shock travels faster than the 
sound speed in the preshock material) then M s can be calculated from Eq. (10): 

M s = M Fe + / dr4vrr 2 H r~ 3 = M Fe + 4vr H In (R s /R Fe ) . (17) 

J RFe 

Rp e is the radius of the stellar iron core before collapse, Mp e the iron core mass, and R s 
the shock position as given by Eq. (14). Combining Eqs. (16) and (17) and plugging in 
Eqs. (11), (14), and (15) our final result is 

2GM 




[£ sn / 4vr# ] 3/2 [t +r cut /c s ] 

, X N (18) 

M Fe - M + 4tt H In ( v s ^ [t + r cut /c 

AM acc decreases with time, therefore early collapse to a black hole or rapidly decreasing 
neutrino fluxes and neutrino heating (i.e., small to) favor massive fallback. Since the 
matter behind the shock is radiation dominated (see Bethe 1993b for a discussion) and 
therefore the adiabatic index is 7 = 4/3, one can easily estimate the sound speed c s 
by again using the assumptions of constant pressure, Eq. (9), and constant density, 
P s ^(M s -M)/(4ttR 3 J3): 

pXl/2 / 4 £ N l/2 



-<7J K Uw-M) ' • (19) 
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Typical masses (M s — M) are 0.2-1 Mq for the times of concern, yielding sound veloci- 
ties of about 0.5-1 X 10 9 cm/s, which is only slightly less than the shock speed (Eq. (11)). 
In order to derive an upper limit of the accreted mass we therefore use the lower bound 
for the time in the denominator of Eq. (18), to + r cut /c s > to, and replace the sound 
travel time r cut /c s in the numerator by to which certainly represents an optimistic case 
because ln(x)/x becomes maximum for 2.72. Plugging representative numbers 

into Eq. (18), M = 1.5 Mq, E sn = 10 51 erg, and Rp e = 3000 km, and taking a favor- 
able, large value of the difference of stellar iron core mass and initial neutron star mass, 
(Mp e — M) ~ 0.2 Mq, one ends up with a maximum mass accreted within a period of 
seconds of 

AM acc < 0.16 M Q (20) 

when the fallback is initiated at a time to = 2 s. 

Although the inequalities in Eq. (8) and Eq. (20) were derived on grounds of very 
optimistic assumptions, the fallback mass needed to account for the three late neutrino 
events in Kamiokande seems hard to reach for the typical conditions in the collapsed 
stellar core of an 18 Mq star. Our estimate of AM acc depends on the quantity H which 
varies with the mass of the progenitor star. For more massive stars H is larger. Yet, the 
estimate of Eq. (18) was obtained by neglecting pressure gradients and, in particular, 
by assuming constant density inside the shock radius R s . This certainly overestimates 
the mass close to the newly formed neutron star and thus the fallback. Numerical 
simulations show that the ejected mass is moving in a thick, dense shell behind the 
outward propagating shock while the neutron star is surrounded by a low-density, high- 
entropy region ("hot bubble"). Accordingly, the fallback turns out to be more than an 
order of magnitude smaller than the limit of Eq. (20) in numerical models that include 
the formation of the hot bubble and the neutrino wind from the protoneutron star (Keil 
and Janka, unpublished). Interestingly, the mass accretion rate in these models peaks 
at about 10 s after core bounce and shock formation. 

It must be emphasized once again that both the lower bound on the required ac- 
cretion mass, Eq. (8), and the upper limit on the possibly accreted mass, Eq. (20), were 
obtained by using extreme assumptions throughout. The numerical value in Eq. (8) 
assumes that gravitational redshift is negligible, that the conversion efficiency of gravi- 
tational binding energy into neutrinos is 100%, and that only v e and v e are emitted and 
carry away equal amounts of energy. Equation (20) on the other hand is too optimistic 
because to = 2 s is the absolute minimum for the time when the re-collapse of some 
matter might have been initiated in SN 1987A, and, moreover, because the density and 
pressure gradients between neutron star and supernova shock and the outward push by 
the neutrino wind were disregarded. Not even for all these extreme assumptions do the 
results of Eqs. (8) and (20) get close. 

However, the analytical considerations in this paper and the numerical studies by 
Keil and Janka are based on very simplified models. They do so far neither include effects 
due to a reverse shock (which can already form at about 0.5 s after bounce, see Janka and 
Miiller 1995, and might increase the fallback on a time scale of several seconds) nor do 
they take into account multi-dimensional effects like rotation or turbulence. Therefore 
the question of fallback masses and time scales requires further investigations. 
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5. Discussion and conclusions 



Two exemplary studies were reported where the neutrino signal from SN 1987A was used 
to constrain properties of the supranuclear EOS and of neutrino-matter interactions. In 
particular, the duration of the neutrino pulse in the Kamiokande experiment marks a 
minimum duration of the neutrino emission which allows one to exclude scenarios which 
predict a more rapid cooling of the nascent neutron star. Such conclusions rely on the 
correct interpretation of the late 3 events in Kamiokande as a part of the protoneutron 
star cooling emission. Moreover, considering and normalizing changes due to the intro- 
duction of new physics relative to the results of a "standard cooling model" requires the 
acceptance of the standard as a meaningful model of SN 1987A. The standard neutrino 
physics in our models is based on the description of the interaction targets as particles 
of nondegenerate or degenerate, ideal Fermigases at thermal equilibrium, an approach 
which is generally used in supernova modelling. Phase space blocking effects are in- 
cluded in the formulation of the neutrino-matter interactions (Burrows and Lattimer 
1986, Keil and Janka 1995). 

Our standard treatment yields very good agreement with neutron star models cal- 
culated by other groups (Burrows and Lattimer 1986, Wilson and Mayle 1989, Suzuki 
1989) although the numerical description employed by those groups partly differed sig- 
nificantly from the one chosen in our models. Moreover, the predicted neutrino signals 
from models with standard physics are in satisfactory agreement with the neutrino pulse 
from SN 1987A. This is quite astonishing, considered the facts that the neutrino-matter 
interactions in dense matter are poorly understood and there are a number of effects 
that should reduce the cross sections significantly compared to their naive values (e.g., 
particle correlations and screening; spin fluctuations as discussed in this work; smaller 
coupling constants in the bulk nuclear medium and relativistic corrections both for 
nucleon recoil and form factors; Seckel, personal communication). Moreover, the pro- 
toneutron star cooling models are spherically symmetrical, thus do not regard possible 
convective processes which might also accelerate the neutrino loss. The long duration 
of the neutrino emission from SN 1987A indicates that the protoneutron star is rather 
opaque to neutrinos but the current knowledge of the source of this opacity is very 
unsatisfactory. 

In detail our results were the following: 
(i) Formation of hyperons and the associated softening of the supranuclear EOS can 
lead to a delayed collapse of the newly formed neutron star to a black hole if the 
neutron star mass is above the maximum stable mass of the cold, deleptonized EOS. 
Since the instability is primarily triggered by the progressing deleptonization of the 
star, the collapse occurs on the neutrino diffusion time scale which is of the order of 
several seconds. Since we found the lepton-loss time to be shorter than the signal 
duration at Kamiokande we concluded that either the neutron star in SN 1987A 
was below the mass limit of the employed EOS, 1.58 Mq (gravitational mass, or 
1.77 Mq baryonic mass), or that the suggested phase transition and softening of 
the EOS due to the creation of hyperons does not occur. Last but not least we 
found the onset of gravitational instability to happen "quietly" i.e. without a late, 
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luminous outburst of neutrino radiation. 

(ii) The duration of the neutrino pulse from SN 1987A and, to a minor degree, also the 
average event energies and event numbers in the Kamiokande and 1MB detectors 
allowed us also to constrain the magnitude of a possible suppression of the axial- 
vector interactions of neutrinos with nucleons. This suppression was conjectured 
by Raffelt and Seckel (1994) as a consequence of rapid spin fluctuations of the 
nucleons which are caused by frequent collisions in a dense medium and lead to 
an averaging effect and therefore to a reduction of the effective target spin "seen" 
by the neutrino. In contrast to naive estimates based on pertubation theory which 
suggest spin fluctuation rates of the order of 20-50 T for typical conditions in a 
newly formed neutron star, we deduced from our cooling models that the spin 
fluctuation rate should be of order T only. This corresponds to a moderate, global 
reduction of the neutrino opacity by roughly 50%. It also implies that the axion 
emission rate from a collapsed stellar core is considerably lower than previously 
used values. In particular, the bounds on the axion mass derived from SN 1987A 
are weakened by roughly a factor 10 to m a < 10 -2 eV. 

(iii) Since our conclusions rely on the duration of the neutrino emission from the pro- 
toneutron star as marked by the last 3 events at Kamiokande, one has to test 
possible alternative scenarios for their ability to explain these late counts and the 
preceding pause of several seconds in the neutrino detection. An obviously attrac- 
tive scenario could be a fast decline of the neutrino luminosities from the protoneu- 
tron star in a time of 2-5 s and a subsequent "quiet" phase of low neutrino fluxes, 
followed by accretion of fallback material which possibly yields sufficient neutrino 
emission to cause the last bunch of Kamiokande events. The rapid cooling of the 
protoneutron star could be a consequence of neutrino opacities that are lower than 
their standard values or might be caused by convectively enhanced neutrino trans- 
port. Due to the low energies of the three late neutrinos they correspond to an 
appreciable neutrino flux and thus a sizable fallback mass is required. Analytical 
as well as numerical studies show, however, that even with very optimistic assump- 
tions it seems very unlikely that in case of an 18 Mq progenitor star of SN 1987A 
fallback occurs which is massive enough to match the considered scenario. 
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